Ni+^'^Ni fusion reaction calculated with the density-constrained time-dependent 

Hartree-Fock formalism 



O 

o 

(N 
D 



> 

ON 

cn 

o 
o 



A.S. Umar and V.E. Oberacker 

Department of Physics and Astronomy, Vanderbilt University, Nashville, Tennessee 37235, USA 

(Dated: February 2, 2008) 

We study fusion reactions of the ^^Ni+^^Ni system using the density-constrained time-dependent 
Hartree-Fock (TDHF) formalism. In this formalism the fusion barriers are directly obtained from 
TDHF dynamics. In addition, we incorporate the entrance channel alignments of the slightly de- 
formed (oblate) ''^Ni nuclei due to dynamical Coulomb excitation. We show that alignment leads 
to a fusion barrier distribution and alters the naive picture for defining which energies are actually 
sub-barrier. We also show that core polarization effects could play a significant role in fusion cross 
section calculations. 

PACS numbers: 21.60.-n,21.60. Jz 



I. INTRODUCTION 

Radioactive ion-beam facilities enable us to study the 
structure and reactions of exotic nuclei, in particular the 
physics properties of the "terra incognita" of neutron- 
rich isotopes One important aspect of these studies 
is a detailed investigation of the heavy-ion fusion process 
of exotic nuclei. This is crucial not only for superheavy 
element formation, but it will also lead to a better under- 
standing of the effective N — N interactions in neutron- 
rich nuclei and of enhanced correlations present in these 
many-body systems. 

Recently, fusion evaporation cross sections for the 
^^Ni-|-^^Ni system have been measured down to the 10 nb 
level This experiment confirmed and improved the 
earlier data Q for the same system and it extended the 
data to extreme sub-barrier energies, thus providing a 
challenge for the theoretical understanding of the fusion 
process between two open-shell nuclei. The primary ob- 
servation was a hindrance of fusion in the ^"^Ni-l-^^Ni 
system at extreme sub-barrier energies in comparison 
to reactions involving other Nickel isotopes such as the 
^*Ni-|-^*Ni system. Earlier coupled-channels calcula- 
tions 0, Q failed to reproduce the data at the extreme 
sub-barrier energies. 

Various hypotheses were developed for the explanation 
of the fusion hindrance phenomenon. In Ref. [6|] the hin- 
drance was attributed to the differing stiffness of Nickel 
isotopes due to nuclear structure effects. An excellent 
coupled-channels fit to the data was obtained by supple- 
menting the effective N — N force used in the double- 
folding potential with a repulsive core to account for the 
nuclear incompressibility effects at the nuclear overlap, 
thus leading to a shallow potential pocket. On the other 
hand, Refs. 0, [1] suggest that at such low energies the 
inner turning point of the heavy-ion potential is smaller 
than the touching point rt = Ri + R2. Thus the valid- 
ity of the frozen-density approximation used in Ref. [G] 
becomes questionable. These authors have proposed a 
two-step model for fusion in which the effects of neck 
formation are approximately included [3| • 

The theoretical analysis of the fusion data generally in- 



volves the determination of a phenomenological ion-ion 
potential such as the Bass model [l^l , the proximity 
potential [ll|, [H, [Hj [llLor potentials obtained via the 
double-folding method fll, [M [O, [Hi- Subsequently, 
the actual fusion cross section is calculated by either us- 
ing barrier penetration mode ls fiol. [TtI JTol l20j . or the 
coupled-channel method [1, i, [IlUlijill . The latter in- 
cludes various excitations of the target and /or projectile 
using the coupled-channel formalism 0, HI], as well as 
the inclusion of neutron transfer, and can be consistently 
applied at energies above and below the barrier [l^] . Ef- 
fectively, the inclusion of each additional excitation leads 
to a modification of the original inert core ion-ion po- 
tential, resulting in a series of effective barriers. One 
common physical assumption used in many of these cal- 
culations is the use of the frozen density or the sudden 
approximation. In this approximation the nuclear densi- 
ties are unchanged during the computation of the ion-ion 
potential as a function of the internuclear distance. Fur- 
thermore, the effects included in channel couplings are 
usually based on the static properties of the participating 
nuclei, which may accurately represent the early stages 
of the collision process, but are expected to change as the 
two ions strongly interact. While these methods provide 
a useful and productive means for quantifying multitudi- 
nous reaction data it is desirable to include dynamical 
effects and make contact with the microscopic theories 
of nuclear structure and reactions. 

Recently, we have developed a microscopic approach 
for calculating heavy-ion interaction potentials which in- 
corporates all of the dynamical entrance channel effects 
included in the time-dependent Hartree-Fock (TDHF) 
description of the collision process These effects 

include the neck formation, particle exchange, internal 
excitations, and deformation effects to all order, as well 
as the effect of nuclear alignment for deformed systems. 
The method is based on the TDHF evolution of the nu- 
clear system coupled with density-constrained Hartree- 
Fock calculations to obtain the ion-ion interaction poten- 
tial. Preliminary calculations for the ^'^Ni-f^'^^Sn system 
highlighted the importance of dynamical deformation ef- 
fects [251 . [26j . Here we give a completed study of fusion 
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cross sections using this formalism. 

In the next section we will summarize some theoretical 
aspects of the density-constrained TDHF theory along 
with methods to incorporate dynamical alignment into 
our calculations, as well as the method used to calculate 
cross sections from the resulting barriers. In Section IIIII 
we present interesting aspects of the reaction dynamics 
and compare our results with experiment and other cal- 
culations. In section Hv] we summarize our conclusions. 



II. THEORETICAL METHODS 
A. Density-constrained TDHF method 

For the calculations of dynamical potential barriers 
for the ^^Ni-|-^''Ni system we have used the density- 
constrained TDHF (DC-TDHF) method. Further details 
of the method can be found in Ref . , here we give a 
short overview. 

The density constraint is a novel numerical method 
that was developed in the mid 1980's [13, HI] and was 
used to provide a microscopic description of the for- 
mation of shape resonances in light systems [28^. In 
this approach the TDHF time-evolution takes place with 
no restrictions. At certain times during the evolution 
the instantaneous density is used to perform a static 
Hartree-Fock minimization while holding the total den- 
sity constrained to be the instantaneous TDHF density. 
In essence, this provides us with the TDHF dynamical 
path in relation to the multi-dimensional static energy 
surface of the combined nuclear system. Since we are 
constraining the total density all moments are simulta- 
neously constrained. The numerical procedure for imple- 
menting this constraint and the method for steering the 
solution to pTDHF(r,i) is discussed in Refs. ^27, _28]. The 
convergence property is as good if not better than in the 
traditional constrained Hartree-Fock calculations with a 
constraint on a single collective degree of freedom. 

In Ref. [131 we have shown that the ion-ion interaction 
potential is given by 



ViR)=EuciR)-EA,-EA, 



(1) 



where -Edc is the density-constrained energy at the in- 
stantaneous separation R{t), while E'ai and are the 
binding energies of the two nuclei obtained with the 
same effective interaction. We would like to emphasize 
again that this procedure does not affect the TDHF time- 
evolution and contains no free parameters or normaliza- 
tion. In practice, TDHF runs are initialized with ener- 
gies above the Coulomb barrier and in Ref. 2j| we have 
shown that there is no appreciable energy dependence to 
the barriers obtained via the density-constrained TDHF 
method. The separation coordinate R is the distance 
between the centers of mass of the two nuclei. 



B. Fusion for deformed nuclei 

The heavy-ion interaction potential between two de- 
formed nuclei depends on the distance vector between 
their centers-of-mass, R, and on the relative orientation 
of their intrinsic principal axis systems which may be 
described in terms of three Euler angles (a, (3, 7) per nu- 
cleus, i.e. in the most general case we have 



V = y(R, ai,^i,7i,a2,/32 , 72) 



(2) 



The expression for V can be simplified if the intrinsic 
nuclear density distributions are axially symmetric; in 
this case, the potential does not depend on the Euler an- 
gles 71 , 72 which describe rotations about the symmetry 
axes. If we put, for convenience, the distance vector in z- 
direction, R = Rbz, the potential between two deformed 
axially symmetric nuclei has the structure 



V ^V{R,/3i,f32,Aa) 



(3) 



The heavy-ion interaction potential is calculated with the 
density-constrained TDHF method for a given set of ori- 
entation angles /3i,/32, Aa. Fortunately, test calculations 
using the double- folding method described in Ref. [1^ re- 
veal that the dependence on the Euler angle Aa is negli- 
gible in our case so we put Aa = 0. In Fig. [T] we show the 
definition of the angles (3i and /32 for two oblate nuclei. 





FIG. 1: Shown are the orientation angles /9i and (32 of the 
oblate Ni nuclei with respect to the collision axis. 

In order to calculate the total fusion cross section at 
energy Ec.m. , we first consider the partial cross section for 
given initial orientations (/3i,/32) of the two nuclei which 
is given by 

00 

(j{E,.^.,I3i,P2) = ^Y.^'2L + 1)Tl{E,,^,,I3i,I32) , (4) 



with /cq = •\/2/z-Ec.m.- The fusion barrier penetrabilities 
TLiEc.ra., P2) are obtained by numerical integration 
of the two-body Schrodinger equation using the incoming 
wave boundary condition (IWBC) method [21I [29l| 



-tf d^ 

2/i (ii?2 



L{L + l)h^ 
2fiR^ 



+ V{R,l3i,l32)-E 



where the quantity V{R, Pi, P2) denotes the heavy-ion 
potential obtained via the density-constrained TDHF 
method. For the numerical implementation we have fol- 
lowed the procedure for the coupled-channel code CC- 
FUL described in Ref. [231 . IWBC assumes that once the 
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minimum of the potential is reached fusion will occur. In 
practice, the Schrodinger equation is integrated from the 
potential minimum, i?niin, where only an incoming wave 
is assumed, to a large asymptotic distance, where it is 
matched to incoming and outgoing Coulomb wavefunc- 
tions. The barrier penetration factor, TL^Ecm., Pi, (32) is 
the ratio of the incoming flux at i?min to the incoming 
Coulomb flux at large distance. 
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FIG. 2: (Color online) Dynamic alignment due to Coulomb 
excitation of ^'*Ni. Shown is the orientation probability as a 
function of the Euler angle /9 in a central collision at internu- 
clear distances R = 1500 fm (blue curve) and at 7? = 16 fm 
(red curve). 



the dominant monopole-multipole part of the Coulomb 
interaction, the orientation probability factorizes as fol- 
lows 



sin(/3i)d/3i sin(/32)d/32 sin(/3i)d/3i sin(/32)d/32 

In the special case of no preferential alignment, i.e. all 
initial orientation angles are equally likely, this factor 
reduces to 



d2p(£;,.„,,/?i,/32) 



sin(/3i)d/3i sin(/32)d/32 



noalign 



In Fig. [2] we show the differential alignment probability as 
a function of the Euler angle /3 used in our calculations. 



III. 



RESULTS 



We have carried out a number of TDHF calcula- 
tions with accompanying density constraint calculations 
to compute V{R, [3i, {32) given by Eq. |T]). A detailed 
description of our new three-dimensional unrestricted 
TDHF code has recently been published in Ref [s^ . 
The code was modified to self-consistently generate ini- 
tial states for ^"'Ni with different orientations. For the 
effective interaction we have primarily used the Skyrme 
SLy5 force (33i |. including all of the time-odd terms. In 
this case the ^^Ni nucleus is essentially oblate having a 
quadrupole moment of -0.45 b. This is also confirmed by 
other calculations Is^, • 



Once the partial fusion cross sections (|4]) for given ori- 
entation angles (/3i,/?2) have been calculated, we have 
to take an average over all initial angular orientations of 
both nuclei 



a{E, 



sin(/3i)d/3i / suv{(32)d(i: 



d^P{E,,^,,f3i,p2) 
sin(/3i)d/3i sm{(32)d[32 



<y{Ec.-m., Pi, 02) , 



where d?P{Ec.m., Pi, P2) represents the alignment prob- 
ability for both deformed nuclei. Details of the dy- 
namic alignment formalism are presented in [s^l- We 
give here a brief summary: For a given incident energy 
Ec.m. we carry out a semiclassical Coulomb excitation 
calculation of the dominant collective levels of the de- 
formed nucleus. The energy levels and i5L-transition 
matrix elements for ^"^Ni are taken from experimental 
data [IH: E2+ = 1.346 MeV, Ei+ = 2.610 MeV and 
M{E2,0+ 2+) = -27.0 e fm^ (oblate deformation). 
The Coulomb excitation calculation starts at very large 
internuclear distances (about 1500 fm) when both nuclei 
may be presumed to be in their respective ground states 
and stops at the ion-ion separation distance R{to) (about 
16 fm). The Coulomb excitation amplitudes determine 
the probability distribution of initial orientations. Using 




FIG. 3: (Color online) Potential barriers, V{R, f^i, l32), ob- 
tained from density-constrained TDHF calculations for the 
64f^i_(_64j^j system. The Euler angles (3i and /32 indicate dif- 
ferent orientations of the deformed ®^Ni nucleus. 



All of our TDHF calculations were done at an initial 
energy of Ec.m. — 98 MeV and separation i?(to) — 16 fm. 
As we have reported in Ref. [IJ the potential barriers ob- 
tained from the density-constrained TDHF method are 
not sensitive to the initial energy (above the barrier). 
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We have tested this again by running a few orientations 
at 112 MeV and did not observe any appreciable differ- 
ence. In Fig. [3] we show the barriers obtained for limiting 
orientations of the ^^Ni nuclei. We have also calculated 
these limiting barriers using other effective interactions, 
SkM* 36] and SLy4 [33|, with essentially no difference. 
While most modern Skyrme parametrizations give essen- 
tially the same potential barrier, using an older Skyrme 
force such as the SIII [37[ interaction results in a higher 
barrier and consequently lower fusion cross sections. This 
issue will be discussed later in the manuscript. 

The physical picture which emerges from the barri- 
ers shown in Fig. [3] is that the total fusion cross section 
strongly depends on the deformation phase space. It also 
shows the fallacy of the often-used statement that a cer- 
tain energy is sub-barrier, which stems from spherical 
systems that can be studied using a single barrier. For 
deformed systems this is dependent on the orientation 
of the nuclei. For the ^^Ni-l-^^'Ni system the only truly 
sub-barrier energies are those below the lowest poten- 
tial barrier corresponding to the (3i = P2 = 90° orienta- 
tion, about Ec.m. = 92 MeV. The fusion cross sections 
corresponding to energies above the lowest barrier will 
be dominated by it since above-barrier cross sections are 
much larger than the below-barrier ones. 
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FIG. 4: (Color online) Total fusion cross section as a function 
of -Ec.m.. Shown are the experimental data (filled red circles) 
and the coupled-channel calculation from Ref. Q] (blue dashed 
curve), and fusion cross sections calculated with the density- 
constrained TDHF method using the SLy5 force (solid black 
line) . 

For the calculation of barrier distributions as a func- 
tion of the orientations angles /3i and P2 we have chosen 
an angular spacing of A/3 — 10°. In principle this requires 
four hundred DC-TDHF calculations, which would be 
very time consuming. However, one can show that some 
of the orientations are equivalent to each other. One ma- 
jor assumption we have made is to assume the equality of 



the angular intervals (0, 7r/2) and (7r/2, tt), which is not 
exactly correct when both nuclei are deformed. In or- 
der to assure that this approximation does not effect the 
lowest energy cross sections (primarily determined by the 
lowest barrier) we have explicitly calculated those angles 
that would have appreciable contribution at these ener- 
gies. In total we have computed twenty potential barriers 
corresponding to various orientations for the SLy5 force. 
In principle even this may not be necessary since all of 
the barriers must fall between the limiting cases shown 
in Fig. [21 Although actual calculations show that a con- 
stant angular interval A/3 does not always lead to equally 
spaced barriers such an extrapolation has a minimal ef- 
fect on the actual cross section calculations. We have 
confirmed this by generating such barriers from the lim- 
iting barriers using a numerical averaging procedure and 
calculating the fusion cross section. This method was 
used for the calculation of the cross sections for the SIII 
interaction. 
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FIG. 5: (Color online) Total fusion cross section as a function 
of -Ec.m.. Shown are the experimental data (filled red circles), 
density-constrained TDHF cross sections using the SIII force 
(dashed black line) , and the density-constrained TDHF cross 
sections using the core orientation with the SLy5 force for the 
lowest energy cross sections. 



In Fig. [3] we show the total DC-TDHF fusion cross 
section as a function of the center-of-mass energy (solid 
black curve) using the SLy5 force. Also, shown are 
the experimental data (filled circles), and the coupled- 
channels calculations of Ref. 2'] (dashed blue line) . Re- 
sults for SLy4 and SkM* interactions are indistinguish- 
able from the SLy5 result. We observe that the DC- 
TDHF calculations, which contain no parameters or nor- 
malization, accurately reproduce the fusion cross sections 
for all energies except for energies -Ec.m. < 87 MeV. We 
believe that the small deviations at other energies are 
largely due to the symmetry assumptions made in align- 
ment averaging. As we stated earlier this is not the case 
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FIG. 6: (Color online) Shown are the density contours in the x — z plane for the TDHF time-evolution of the ®'*Ni-|-^'*Ni 
system. Initially both nuclei are oriented with angles f5i = {32 = 90° at a center-of-mass energy of Ec.m. ~ 98 MeV. The 
values for the ion-ion separation R correspond to special points along the potential barrier at the lowest experimental energy 
of -Ec.m. ~ 86 MeV; (a) outer turning point, (b) inner turning point, (c) reorientation of the core, and (d) potential minimum. 



for the lowest energies as these were explicitly done with- 
out symmetry assumptions. 

In Ref. @ the hindrance of the fusion cross section 
at deep sub-barrier energies was attributed to nuclear 
incompressibility effects at the nuclear overlap and in 
Ref. |8)] to the modifications of the inner turning point 
in comparison to standard folding model calculations us- 
ing the frozen-density approximation. As we shall see 
below, both of these assertions may be correct. In or- 
der to investigate the dependence on incompressibility 
we have repeated our calculations using an older Skyrme 
parametrization, SIII (37i | . having a nuclear matter in- 
compressibility coefficient of 356 MeV, as opposed to 
most of the modern Skyrme parametrizations that have 
nuclear matter incompressibility around 230 MeV. Al- 
though the SIII interaction accurately reproduces many 
spectroscopic properties of spherical and close to spher- 
ical nuclei it is known not to correctly reproduce fission 
barriers and other phenomena involving large deforma- 
tions. In Fig. [5] the black dashed line corresponds to the 
DC-TDHF calculations using the SIII interaction. As we 
see, this interaction underestimates the fusion cross sec- 
tions at intermediate energies but does relatively well at 
low and high energies. Although it is difficult to disen- 
tangle the effect of incompressibility from all the other as- 
pects of the effective interaction, one may conclude that a 
new Skyrme force with a higher than the accepted nuclear 
matter incompressibility of 230 MeV and fitted to repro- 



duce large amplitude collective phenomena could better 
reproduce the fusion cross sections. 

In order to better examine the evolution of the nu- 
clear density, in Fig. [S] we have plotted the nuclear 
density at four special internuclear distances R for the 
Pi = P2 = 90° initial orientation. Frame (a) corresponds 
to the nuclear density at the outer turning point of the 
ion-ion potential (i? = 13.1 fm) at Ec.m. < 86 MeV. One 
striking observation from this frame is that the orien- 
tation of the nuclear core seems to be rotated in com- 
parison to the total nuclear density by 7r/2. Frame (b) of 
Fig. [S] shows the total density at the ion-ion separation of 
R — 1 1.0 fm, which approximately corresponds to the lo- 
cation of the inner turning point at Ec.m. < 86 MeV. The 
orientation of the nuclear core is still the same as at the 
outer turning point. At around i? = 10 fm the nuclear 
core rotates and aligns with the total nuclear density. In 
frame (c) we show this at i? = 9.5 fm. The last frame 
(d) shows the nuclear density at the potential minimum 
which occurs around R = 8.0 fm. Based on the above ob- 
servation we can make a conjecture that while at higher 
energies the potential barrier is largely determined by 
the nuclear surface at deep barrier energies core nucle- 
ons play a significant role in barrier dynamics. This may 
manifest itself as repulsive core in constructing a poten- 
tial model for the problem. In our calculations we can 
incorporate this effect by using the alignment angles of 
the core rather than the total density in the calculation 
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of the cross sections at the lowest energies. This is shown 
by the black solid curve in Fig. O As we see the results 
are in excellent agreement with the data at low energies. 

IV. CONCLUSIONS 

As we investigate fusion reactions involving neutron 
rich and deformed nuclei it is apparent that an under- 
standing of the structure of these nuclei is crucial to 
the description of the reaction dynamics. For these nu- 
clei various effects, such as inelastic excitations, particle 
transfer, and other dynamical effects lead to substantial 
modification of the naive potential barrier calculations 
which assume an inert core and no dynamics. Conse- 
quently, the definition of sub-barrier fusion becomes am- 
biguous since it is difficult to determine the barrier a 
priori. 

We have performed density-constrained TDHF calcu- 
lations of fusion cross sections for the ^'*Ni-|-^''Ni system. 
Our results agree well with the measured data despite 
having no adjustable parameters. This indicates that 
many of the reaction dynamics are included in the TDHF 
evolution of the nuclear density. We have also investi- 
gated the dependence of our results on the microscopic 
effective interaction. We find that while all of the mod- 
ern Skyrme parametrizations show very small deviations 
in the fusion cross sections, which was also observed in a 
more systematic study of spherical systems ^35] , the older 
parametrizations yield very different results. Since older 



Skyrme parametrizations were fitted mostly to reproduce 
properties of spherical nuclei they may not be a good can- 
didate for fusion studies. On the other hand it may be de- 
sirable to investigate the dependence on incompressibility 
using a more modern Skyrme parametrization, which is 
not yet available. 

We have further investigated fusion cross sections at 
deep sub-barrier energies. In the absence of a true many- 
body tunneling approach to nuclear fusion it is difficult to 
envision the dynamical formation of the potential barrier 
at very deep sub-barrier energies. As we go further down 
in energy the inner turning point of the ion-ion potential 
involves larger overlaps between the participating nuclei. 
Consequently, core nucleons may play a more dominant 
role in dynamically building up the potential barrier. In 
the case of f3i — (32 — 90° corresponding to the low- 
est potential barrier we observed that the nuclear core 
has a different orientation from the total nuclear density. 
Based on this observation we have speculated that at the 
lowest energies it may make sense to use the orientation 
of the core rather than the nuclear surface. This core po- 
larization effect allows us to reproduce the experimental 
cross sections at the lowest energies. 
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